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Abstract 

We construct a three-dimensional analytically soluble model of the nonlinear 
effects in Coulomb excitation of multiphonon Giant Dipole Resonances (GDR) 
based on the SU(2,1) algebra. The full 3-dimensional model predicts further 
enhancement of the Double GDR (DGDR) cross sections at high bombarding 
energies. Enhancement factors for DGDR measured in thirteen different pro- 
cesses with various projectiles and targets at different bombarding energies 
are well reproduced with the same value of the nonlinearity parameter with 
the exception of the anomalous case of 136 Xe which requires a larger value. 
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One of the most interesting applications of Coulomb Excitation in heavy ion collisions 
]1|-[8|] is investigation of multi phonon nuclear Giant Resonances (GR) 0. Possibility to ex- 
cite multiple GR involves Bose statistics of collective excitations and constitutes the "family 
property" of vibrational collective motion in both finite (nuclei, clusters) and infinite quan- 
tum systems. Within this concept, making no distinction between infinite and finite systems, 
the excitation process can be modeled via a multidimensional quantum oscillator coupled lin- 
early to the external time-dependent field, providing excellent agreement with the single-GR 
experimental data ||. 

Validity of this completely linear theory has been questioned by nearly all the exper- 
imental data wherever multi-phonon GR has been observed @: the double Giant Dipole 
Resonance (GDR) excitation cross sections || are found about 1.3 — 2 bigger than follows 
from theory. This shortcoming of the linear theory, known as "enhancement factor prob- 
lem" , has been addressed widely in current literature within a number of approaches: higher 



order perturbation theory M, anharmonic effects M, Jul, concept of hot phonons [11-13 



Clearly, nonlinear effects, that in principle can not be neglected in finite Fermi systems [H 
are not easily dealt with either at a microscopical or even at a phenomenological level |7]]. 

It is therefore appealing to construct a natural, soluble "minimal extension" of the har- 
monic model of Coulomb excitations, in which deviations from the linear scheme are rea- 
sonably described via a few parameters. Without dealing in depth with the microscopic 
theory, we present here such a "minimal extension" . This single-parameter nonlinear model 
can be solved exactly using algebraic properties of boson operator combinations forming 
algebra SU(2,1). The model allows us to correlate all the experimental data for the "en- 
hancement factors" in various nuclei and various bombarding energies, using a single value 
of the universal nonlinear parameter. 

Within the semiclassical approach [[HJ to Coulomb excitation, the projectile motion is 
approximated by a constant velocity v on a straight line classical trajectory with impact 
parameter b and internal excitation is treated quantum mechanically. The intrinsic state 
|\I/(t)) of the system undergoing excitation obeys the time-dependent Schrodinger equation, 



id/dt\m{t)) = [H + v(t)} \*(t)), |*(* = -oo)) = |o>, (l) 

where H is the intrinsic Hamiltonian and V(t) = — T>\^\ + Vo(t)T>l + h.c. is the 

channel-coupling interaction with T>^ and T> the dipole operators acting in the space of the 
multi-GDR states created by the boson operators d+ , with the angular momentum projection 
m (h=c= 1). The functions t>i(o)(£) describe the interaction with the electromagnetic field 
fl> 01 ( see below). The excitation probability of an intrinsic state \N) with N phonons in 
a collision with impact parameter b and the total cross section <jn are PJ^,p|, |l5|1 : 

/oo 
bP N (p). (2) 



The internal nuclear Hamiltonian H is nearly harmonic |2[, ||16|| , |17| , so H = uoN = 
ujJ2 m d m d m . In principle, this does not exclude anharmonicities |Tj| in the transition oper- 
ators T>\ T> when expanded in terms of phonon operators 

= d m + x' s ^d\ n d\ ni d mi + '^x mi d\ n d\ ni d\ ni + x 2 E d^d^d^d^d^ + (3) 



These effective nonlinearities [H|, [[| could result from perturbation theory treatment of 
anharmonicities in the phonon Hamiltonian, from coupling to other degrees of freedom, both 
collective (e.g., quadrupole GR) and noncollective JTT| etc. The linear limit of the problem, 



Tr m = d^, is exactly soluble giving the Poisson formula for the excitation probabilities 

P N = e~ pP Wv P= E l«r m ! 2 = E 1/ v m (t)e^dt\ 2 (4) 

iV ' m=0,±l m=0,±l J -°° 

where the amplitudes a^ arm are given by the modified Bessel functions Ki and K [0. 

In order to reduce the number of unknown parameters in Eq.(|3D it is reasonable to 
restrict the higher-order nonlinear corrections ocxi. We save the first, dominating, nonlinear 
term in ([3D with its coefficient CKx-Cl and save only those terms in ([3D which match terms 
appearing in the expansion of the square root, V^ m — > d^ m {l + 2xY, m d m d rn ) 1 /' 2 . 

The nonlinear effects are now controlled by the single parameter x. This may be justified, 
by noting that the leading term is saved, while the ansatz for the higher order terms, which 
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should be small anyway, obeys the basic requirement that they get smaller as the order 
increases. Further, this parametrization leads to a soluble problem. 

The group theoretical solution is based on the consideration of the eight operators 

1 



D + = (d± 1 -d+ 1 ) y /k + N/2, D+ = d£V2k + N, J 

1 r 
4 



(dli - <x)(d-i - d+i) + 2(2k + N) 



2V2 
1 



d_ x — ) d\ 



'0- 



2k + N + dJd 



(5) 



and D~=(D + Y , Z)g — (-^o") ; (J + V with the Casimir invariant. They form a 

closed SU(2,1) algebra, the non-compact analogue of SU(3). The pairwise commutators be- 
tween (|) can be evaluated directly: Three of them are [D~ ,D + ]=2D°, [Dq ,Dq ]=2Dq, 
[J~,J + ]=2(D° -D°). The other nonzero ones are [D~ ,D°]=2[D~ ,D%]= [D ,J-]=D-; 
[D ,D° }=2[D ,D°}= [D-,J+]=D , and [D~ ,D^]=2[J~ ,D°]= -2[J~ ,D%]=J~ . The remain- 
ing nonzero commutators are given by Hermitean conjugates to the above. In the interaction 
representation, the evolution equation i-^\ip(t)) = e lHot V(t)e~ lHot \^(t)) ([!]) with (|) reads: 



i{d/dt)\m) = 2x i/2 \vi(t)& + vt(t)D- + v (t)Di + v* (t)D } \m)- 



(6) 



where v\ = v\e wt and vq = voe lult /V2. Any product of exponentials involving the operators 
from the set (|5|) can be reduced to a simpler exponential (see, e.g., as due to closure of 
the pairwise commutators between fl5|), no new operator structures arise while re-arranging 
order of the operators in the time-ordered exponential (see, e.g., fl20f ) that solves (0),® 

|^(t)> = Texpif e iH ° T V(T)e- iH ° T dr} |0) = e aD++bD ° e ^ +dDl e9 j- e a' D - + b'D-^ (7) 

where the eight time-dependent c-numbers, the Latins a — g, must be chosen so if)(t)) (0) 
obeys the Schrodinger equation @. The ansatz (|7|) can be simplified using D~\0) = Dq\0) = 
J^O) = 0; D°\0) = Dq\0) = k\0), which follow from (|5|). The expression (|7|) then reduces to 



4x (\a(t)\ 2 + P(t)\ [ 



D, 



°|0), 



(8) 



where the first factor comes from unitarity, (ip(t)\ip(t)) = 1. The phase (f)(t) is unimportant 
for the following. Substituting (H) into (fj|), we obtain, after some lengthy algebra using 



heavily the operator identity e A Be A = B+ [A, B] + (1/2!) [A, [A, B]] + the two nonlinear 
equations for the amplitudes a and (5 



ida/dt — Qa = i'i(t), id/3/dt — Q(3 = v (t), 



(9) 



with Q = 4x[-Ug(t)/3 + v*(t)a] and the initial condition a(— oo) = /?(— oo) = 0. After 
projecting the asymptotic state (§) at t = oo onto the states with definite number of GDR 
phonons, we obtain the non-Poissonian expression for the excitation probabilities 



N 



T(± + N) [4x(\a\ 2 + \ffl 



N 



(10) 



N\r{±) [l-4x(|a| 2 +^P 

where a and (3 are the asymptotic solutions of the system (§) at t = oo. Their values 
can be easily tabulated by solving (|j). If x — > 0, Eq . (p^0|) is reduced to Poisson, while the 



solutions to (H) are reduced to the harmonic amplitudes (f|), \a\ 



and |^| 



a 



harm 




|/\/2 = ^-i^o(0- Thus, the harmonic results are restored. At nonzero 



nonlinearity x > 0, the multiple GDR excitation probabilities Pjv turn out to be enhanced 
as compared to their values p^ arm in the harmonic limit, Eq.(|4]). The cross sections o^r 
and their harmonic values, a^ rm are given by integrating in Eq.([|) from the grazing value 
1.2(v4g/ 3 + A 1 ^). The enhancement factors r N = o^ja^^ 1 can be studied and compared 



with the data. The functions v m {t\ in the notations £ = — , r = are 



vi(t) = Ff 



v 



(t) = FV2 



7 



df 



- 



Z sp e 2 j 
2b 2 



A 2 J x 3 m N ■ 80MeV 



Or 

where / = (1 + r 2 )~^ and 7 = (1 — v 2 )~^. In the strength F, ttin and e are the proton 
mass and charge, Z, N and A denote the nuclear charge, the neutron number and the mass 
numbers, the labels ex (sp) refer to the excited (spectator) nucleus in the colliding pair. 

Let us discuss the energy dependence of the two-phonon enhancement factor r 2 . In Ref. 
PH , we have studied this behavior for the Pb+Pb system at the bombarding energies in range 
70 — 700MeV-A, using truncated (two-dimensional) dynamics and neglecting the longitudinal 
response in the excitation process (vq = 0). Within this "toy model" |2TJ, the dynamics is 
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described by the three operators D,D + and D° which form the SU(1,1) subalgebra of the 
SU(2,1). The enhancement factor drops steadily as the bombarding energy grows. This is 
not the case at higher energies, when the transverse approximation breaks down and full 
solution based on the SU(2,1) algebra is required. This is illustrated in Fig. la, where the 
enhancement factor r 2 is shown as a function of 7 for the case of Pb + Pb collision in the 
energy range up to 4 Mev, as compared to the "transverse approximation" . 

The remarkable fact about the full 3-dimensional model is that the enhancement factor 
tn starts to grow again as the relativistic factor 7 passes the "extremal" point (~ 1.3 — 1.5 for 
heavy colliding nuclei). Thus the 3-dimensional SU(2,1) model predicts a new interesting 
qualitative effect which can be tested in experiments using higher bombarding energies. 
This behavior of r N is related to the 7-dependence of the "scaled" longitudinal function v 
(Eq.|TTD v {t) = v {t)e iu4 = F^2 [-7^ (/e^ T ) +i^fe^ T ]. The 7-dependence of the 

two terms here is very different. While the first term scales as 7, the absolute maximum of 
the second one scales as £/7 = r y~ 2 . The first term is proportional to the time derivative of 
the function ^^2^3/2 which dies out at £=±00. Therefore, in the harmonic limit, when the 
longitudinal amplitude (3 reduces to the simple integral over time of the function v , the first 
term simply vanishes. The harmonic solution is then given by the integral of the second term 
to give oc -ifo(£) which dies out as I/7 as bombarding energy grows, and the longitudinal 
response becomes negligible at high bombarding energies as is well known 0-44]. This is 
not the case at nonzero nonlinearity x 7^ 0: the amplitude (3 is given by the solution to the 
coupled nonlinear system (|^), and the first term in v not only does contribute but in fact 
becomes dominant. Unlike the second term, the first one grows with 7 and this determines 
the behavior of the enhancement factor. 

We present below in Figs. lb, Fig.2 and Fig.3a exact numerical results for the cross 
sections calculated according to Eqs.(0), and ([K]) and using numerical solution of Eq.(^). 
Since the nonlinear parameter x is an internal property of the nucleus in which the GDR is 
excited, it is reasonable to expect that it varies from one nuclear species to another. For the 
sake of comparison with the experimental data, it is expedient however to choose a single 
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average value x = 0.29 for all cases and use it to calculate the cross sections. As seen in 
Figs.la,2 and Fig. 3a, the nonlinear model can in this way (no fitting) reproduce rather well 
the experimental values || of r2 for the twelve excitation processes of Pb (Fig.2a,b) and Au 
(Fig. 2c), using different bombarding energies and targets/projectiles. An exceptional case is 
Xe, where considerably larger value of x is required. To illustrate this, the optimal values of 
x for each individual datum are shown in Fig. 3b, together with the adopted average value. 

In conclusion, we presented a simple soluble model to account for the nonlinear effects 
in the transition operators for the Coulomb excitation of multi-phonon GDR via relativistic 
heavy ion collisions. The solubility of the model is based on the group theoretical properties 
of the boson operators. It allows to construct the solution for the dynamics of the multi- 
phonon excitation within the coupled-channel approach. The well known harmonic phonon 
model appears to be a limiting case of the present model when the nonlinearity goes to zero. 
The main advantages of the limiting harmonic case (unrestricted multiphonon basis, preser- 
vation of unitarity and possibility of analytical treatment) remain present in our nonlinear 
scheme. Therefore, the model can be viewed as a natural extension of the harmonic phonon 
model to include the nonlinear effects in a consistent way while keeping the model solvable. 
With proper modifications, it can be used in other (nuclear, molecule etc.) problems. 

At low enough bombarding energies, the enhancement factor drops as the bombarding 
energy grows. This is consistent with the data and gives results similar to those recently 
obtained in a different context, with a theory based on the concept of fluctuations (damp- 
ing) and the Brink- Axel mechanism |pT| , ||12|| . The interesting property of the full three- 
dimensional results obtained here is that the enhancement factor starts growing again at 
high bombarding energies (~ 0.3— 0.5GeV- A for heavy colliding nuclei.) Besides being an 
interesting prediction to be tested in experiment, this behavior allows one to separate the 
effects of nonlinearity considered here and the effects proposed in works ||iT|-[i"3|j. 

The work has been supported by the FAPESP and by the CNPq. 
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Figure Captions 

Fig.l. (a) The enhancement factor r 2 as a function of 7 for 208 Pb+ 208 Pb collision. The 3D 



results vs. the transverse approximation results (SU(1,1) toy model, for x = 0.19. 

(b) Enhancement factor r 2 = a 2 /a% arrn for the DGDR excitation in 208 Pb SIS db function 
of 7 (open symbols connected by a line). The data || (filled symbols): triangle up - Zn 
projectile, triangle right - U tagret, circles - Pb target and projectile. The value of the 
nonlinear parameter is kept fixed, x = 0.29. 

Fig. 2. (a) The same as in Fig.lb, x = 0.29. Square - Ar projectile, diamond - Kr 
projectile, triangle down - Ho, and circle - Sn. (b) The same for the DGDR in Au. Circle 
- Kr projectile, square - Au projectile, diamond - Bi projectile and triangle - Ne projectile. 
Open circles - theoretical values for the Bi projectile (the results for Kr and Au projectiles 
are the same at 7 ~ 2). 

Fig.3. (a) The same as in Fig.lb,Fig2, but for the DGDR in Xe, x = 0.29. 

(b) The optimal values of the nonlinear parameter x for the thirteen processes vs strength 
parameter F/u, Eq. (|TT|) . The processes are numbered in the order they appear in Figs. lb, 
2a, 2b, 3a (from left to right). 
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